164 


IEEE  GEOSCIENCE  AND  REMOTE  SENSING  LETTERS,  VOL.  3,  NO.  1,  JANUARY  2006 


A  Statistical  Approach  to  WindSat  Ocean 
Surface  Wind  Vector  Retrieval 

Craig  K.  Smith,  Member,  IEEE ,  Michael  Bettenhausen,  Member,  IEEE ,  and  Peter  W.  Gaiser,  Senior  Member,  IEEE 


Abstract — WindSat  is  the  first  space-based  polarimetric  mi¬ 
crowave  radiometer.  It  is  designed  to  evaluate  the  capability  of 
polarimetric  microwave  radiometry  to  measure  ocean  surface 
wind  vectors  from  space.  The  sensor  provides  risk  reduction  for 
the  National  Polar-orbiting  Operational  Environmental  Satellite 
System  Conical  Scanning  Microwave  Imager/Sounder,  which 
is  planned  to  provide  wind  vector  data  operationally  starting 
in  2010.  The  channel  set  also  enables  retrieval  of  sea  surface 
temperature,  and  columnar  water  vapor  and  cloud  liquid  water 
over  the  oceans.  We  describe  statistical  algorithms  for  retrieval  of 
these  parameters,  and  a  combined  statistical/maximum-likelihood 
estimator  algorithm  for  retrieval  of  wind  vectors.  We  present  a 
quantitative  analysis  of  the  initial  wind  vector  retrievals  relative 
to  QuikSCAT  wind  vectors. 

Index  Terms — Microwave  radiometry,  polarimetry,  retrievals, 
wind. 


I.  Introduction 

INDSAT  is  a  conically  scanning  polar-orbiting  multi- 
frequency  microwave  radiometer  launched  in  January 
2003.  Its  objective  is  to  measure  the  partially  polarized  emis¬ 
sion  from  the  ocean  surface  and,  therefore,  test  and  fully  eval¬ 
uate  the  viability  of  using  polarimetric  microwave  radiometry 
to  retrieve  ocean  surface  wind  vectors.  WindSat  measures  po¬ 
larized  radiometric  brightness  temperatures  (Tg)  in  16  “chan¬ 
nels”:  vertically  and  horizontally  polarized  Tg  (Tv  and  Th)  at 
6.8,  10.7,  18.7,  23.8,  and  37.0  GHz,  and  third  and  fourth  Stokes 
parameters  (Tjj  and  T4)  at  10.7, 18.7,  and  37.0  GHz.  Tjj  and  T4 
are  derived  from  the  difference  of  ±45°  linear  and  left/right  cir¬ 
cular  polarizations,  respectively.  The  WindSat  reflector  antenna 
has  an  1 1-feedhom  cluster,  resulting  in  fixed  nominal  Earth  in¬ 
cidence  angles  (EIA)  ranging  from  50°  to  56°  (all  horns  within  a 
frequency  band  have  the  same  EIA)  [1].  The  El  As  are  known  to 
better  than  0.05°  [2].  A  detailed  description  of  the  radiometer 
can  be  found  in  [1].  WindSat  provides  risk  reduction  for  the 
Conically  Scanning  Microwave  Imager/Sounder  to  be  flown  on 
the  National  Polar-orbiting  Operational  Environmental  Satellite 
System  [3]. 

We  have  developed  an  empirical  statistical  regression  tech¬ 
nique  for  retrieval  of  sea  surface  temperature  (SST),  wind  speed 
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TABLE  I 

Nominal  NEDT  Values  for  Rain-Free  Ocean 
Retrieval  Cells  (in  Kelvin) 


6.8 

10.7 

18.7 

23.8 

37.0 

H 

< 

H 

0.111 

0.036 

0.033 

0.063 

0.023 

TuT4 

- 

0.053 

0.049 

- 

0.036 

(W)  referenced  to  10-m  height,  compass  wind  direction  (cp), 
columnar  water  vapor  (V),  and  cloud  liquid  water  (L).  The  pur¬ 
pose  of  this  effort  is  threefold:  1)  obtain  a  quick  initial  assess¬ 
ment  of  WindSat’ s  wind  vector  measurement  capabilities;  2) 
provide  a  retrieval-oriented  test  platform  for  enhancements  to 
the  WindSat  Temperature  and  Sensor  Data  Record  algorithms; 
and  3)  provide  either  an  independent  set  of  retrievals  from,  or 
a  priori  data  for,  a  physically  based  algorithm.  The  regression 
technique  and  the  selection  of  regression  training  data  were  de¬ 
signed  to  facilitate  retrieval  of  all  five  parameters,  as  discussed 
in  Sections  II  and  V.  However,  in  this  letter  we  focus  on  appli¬ 
cation  of  the  regressions  to  wind  vector  retrieval. 

II.  Regression  Technique 

For  this  initial  validation  of  the  regression  technique,  we  used 
WindSat  Tg  averaged  to  the  lowest  resolution  (6.8  GHz),  where 
the  Tg  for  all  16  channels  are  available.  The  Tg  for  each  channel 
are  sampled  at  approximately  12.5  km  along  scan  and  along 
track  and  averaged  to  a  common  spatial  resolution  of  approx¬ 
imately  40  x  60  km.  The  nominal  measurement  noise  for  the 
retrieval  cells,  specified  by  the  noise-equivalent  differential  tem¬ 
perature  (NEDT)  for  mean  ocean  scene  Tg,  is  given  in  Table  I. 

The  regression  training  (i.e.,  regression  coefficient  determi¬ 
nation)  is  performed  empirically,  using  collocations  of  WindSat 
Tg  with  various  sources  of  satellite  retrievals  and  numerical 
weather  prediction  model  (NWP)  data  (Section  V).  Unlike  re¬ 
gressions  trained  using  Tg  simulated  with  a  physical  model  [4], 
[5],  empirical  regressions  automatically  compensate  for  bias  er¬ 
rors  in  the  measurements,  and  potential  problems  such  as  time- 
dependent  calibration  errors  or  high  noise  results  in  deemphasis 
(coefficients  near  zero)  for  the  affected  channel.  Therefore,  we 
use  all  16  channels  in  all  regressions.  Regression  coefficients  are 
determined  by  linear  least  squares. 

We  use  a  linear  method  to  retrieve  parameters  for  which  the 
relationship  with  the  Tg  ranges  from  nearly  linear  to  severely 
nonlinear.  Three  techniques  exist  to  address  this  nonlinearity 
problem.  The  first  uses  functions  of  the  Tg  that  are  linear  with 
respect  to  a  given  retrieval  parameter  as  arguments  in  the  regres¬ 
sions  [4],  [5].  Additionally,  terms  that  are  quadratic  or  higher 
order  in  the  Tg  or  the  linearization  functions  may  be  included 
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to  allow  the  regression  training  to  eliminate  some  of  the  nonlin¬ 
earity  directly  [5].  The  second  technique  is  to  retrieve  a  func¬ 
tion  of  the  desired  parameter,  on  which  the  Tg  have  an  approx¬ 
imately  linear  dependence  [5].  The  third  method  is  to  train  and 
apply  regressions  on  restricted  domains  (“bins”)  of  a  geophys¬ 
ical  parameter  on  which  the  Tg  have  a  nonlinear  dependence. 
The  bins  are  chosen  so  that  the  dependence  of  the  Tg  on  that 
parameter  is  approximately  linear  within  each  bin  [4] . 

The  application  of  these  techniques  to  the  retrievals  is  ex¬ 
plained  below.  We  make  a  significant  extension  to  the  second 
technique  to  retrieve  wind  direction,  for  which  the  signal  in  all 
the  Tg  is  severely  nonlinear  (in  fact  nonunique)  and  depends  on 
wind  speed.  There,  we  perform  regressions  for  several  functions 
of  wind  speed  and  direction,  and  then  use  the  regression  results 
in  a  maximum-likelihood  estimator  (MLE)  for  wind  direction 
retrieval,  as  discussed  in  Section  III. 

The  regression  for  any  retrieval  parameter  P  is  of  the  form 

16  16 

P  =  CLo  +  ^  ^  CLiti  +  ^  ^  bitf  +  Cie  +  C2C2  (1) 

i=  1  i—i 

where  p  is  the  WindSat  brightness  temperature  for  channel 
i(Tgi),  except  for  the  23.8-GHz  channels,  where  ti  = 
ln(290  —  Tgi);  the  latter  serves  to  linearize  U  with  respect  to 
the  V  signal  when  water  vapor  absorption  is  appreciable  [4], 
[5],  which  significantly  improves  the  V  retrieval  but  does  not 
significantly  affect  any  of  the  other  retrievals.  The  quadratic 
terms  greatly  reduce  the  nonlinear  interactions  (“crosstalk”) 
between  SST,  W,  V,  and  L.  Inclusion  of  the  EIA  for  the  37-GHz 
vertical-horizontal  (vh)  feedhorn  (e)  reduces  the  SST  retrieval 
error  by  5%,  but  has  little  effect  on  the  other  retrievals.  The 
small  variations  from  the  nominal  EIAs  for  the  different  feed- 
horns  show  a  linear  relationship  to  each  other,  so  that  inclusion 
of  a  single  EIA  is  sufficient  for  a  linear  regression.  The  regres¬ 
sions  allow  for  direct  retrieval  of  SST,  W,  V,  and  L. 

The  directional  component  of  the  microwave  sea  surface 
emissivity  (e^ir)  can  be  expressed  as  a  two-term  Fourier  series 
in  relative  wind  direction  ((/?r,  the  compass  wind  direction 
minus  radiometer  look  direction  projected  on  the  Earth)  [6]. 
The  Fourier  coefficients  are  nonlinear  functions  of  wind  speed 
that  depend  on  Earth  incidence  angle,  frequency,  and  polariza¬ 
tion  [6] .  The  series  is  even  for  Tv  and  Th  and  odd  for  Tjj  and 
T4  [6],  [7] 

£dir  =  fi(w)  cos ipR  +  f2(W)  cos 2ipR;(Tv,Th) 

£dir  =  gi{W)smipR  +  g2{W)sm2tpR\(Tu,TA).  (2) 

Thus,  for  ip  retrieval,  we  perform  four  regressions,  for  the 
along-look  wind  speed  (C/i  =  VEcos(/2r),  cross-look  wind 
speed  (U2  >=  IE  sin  <^r),  and  “second  harmonics”  of  the 
along-  and  cross-look  wind  speeds  (C/3  =  W  cos  2<^r  and 
C/4  =  W  sin  2(£>r).  The  along-look  wind  speed  is  identical  to 
the  “line-of-sight”  wind  speed  retrieved  by  Wentz  [8].  While  the 
U\  and  C/2  regressions  are  sufficient  for  unique  wind  direction 
retrieval,  we  find  significant  improvement  for  W  >  9  m/s  by 
using  all  four  regressions  as  input  to  an  MLE  for  wind  direction 
retrieval. 


The  regressions  are  accomplished  using  a  two-stage  tech¬ 
nique.  In  the  first  stage,  the  retrievals  are  performed  for  W, 
V,  and  L,  using  regression  coefficients  trained  on  the  entire 
training  dataset  (Section  V).  In  the  second  stage,  the  first-stage 
W  regression  result  is  used  to  determine  a  wind  speed  bin,  from 
which  regression  coefficients  are  selected  to  perform  the  SST 
and  U\  through  C/4  retrievals.  The  bins  are  2  m/s  wide,  and  cover 
the  wind  speed  ranges  0-2,  2-4,  4-6  m/s,  etc. 

The  second-stage  regressions  offer  significant  advantages 
over  analogous  first- stage  regressions.  Within  each  bin  the 
Fourier  coefficients  in  (2)  are  approximately  linear  functions  of 
W — i.e.,  the  wind  direction  signals  are  approximately  linear  in 
U\  through  C/4 — improving  the  U\  through  C/4  linear  regression 
and  MLE  wind  direction  performance. 

Second-stage  performance  is  improved  by  widening  the  bins 
during  regression  training,  to  account  for  the  “misbinning” 
resulting  from  retrieval  error  in  the  first-stage  W  regression 
[0.8  m/s  root  mean  square  (rms)  error].  However,  the  benefit 
of  accommodating  misbinning  is  counterbalanced  by  the  in¬ 
creasing  nonlinearity  of  the  Tg  with  respect  to  U\  through  C/4 
as  the  bins  are  widened.  Given  the  bin  width  of  2  m/s  for  the 
retrievals,  a  bin  width  of  4  m/s  for  the  training  is  optimal;  the 
bins  are  thus  widened  by  1.0  m/s  on  each  side  for  second-stage 
training  (0-3,  1-5,  3-7  m/s,  etc.). 

III.  Maximum-Likelihood  Estimator 

Following  application  of  the  second- stage  regressions,  U\ 
through  C/4  are  provided  to  the  MLE  for  wind  direction  retrieval. 
The  likelihood  function  to  be  maximized  is  for  estimating  two 
parameters  (W,  <^r)  from  four  “measurements”  (Z7i  through 
C/4)  whose  errors  are  correlated  and  approximately  Gaussian 
distributed  [9].  Therefore,  the  MLE  finds  multiple  wind  vector 
solutions,  (W,pr)  pairs  known  as  “ambiguities,”  as  local 
minima  of 

X2  =  suT  •  c^1  •  su 

where 

/  Ui-W  cos  cpR  \ 
cfr  U2-W  sin  cpR 
°  I  C/3  —  VE  cos  2(pR  * 

\U4-W  sin  2  ipR  / 

Cu  is  the  error  covariance  matrix  of  the  U\  through  C/4  re¬ 
trievals  for  the  bin  determined  by  the  first-stage  W  regression 
(Section  II)  and  is  obtained  by  running  the  regressions  on  the 
training  dataset  (Section  V).  The  SU  vector  is  simply  the  differ¬ 
ence  between  the  U\  through  C/4  regressions  and  their  mathe¬ 
matical  definitions  given  in  Section  II. 

The  MLE  may  appear  to  be  a  two-dimensional  search  for 
chi-square  minima  in  the  space  of  (IE,  <^r).  However,  it  should 
be  noted  that,  because  U\  through  C/4  are  linear  in  W,  mini¬ 
mization  along  the  W  axis  yields  a  single,  analytical  solution 
for  W  at  each  value  of  (/?r,  denoted  VEmin(^R).  Substitution 
of  Wmin((pR)  into  the  definition  of  the  chi-square  results  in  a 
one-dimensional  minimization  along  the  ipR  axis.  This  one-di¬ 
mensional  minimization  is  subject  to  the  following  conditions: 
1)  VEmin(^R)  is  nonnegative  (physical  solution);  2)  the  second 
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partial  derivative  of  the  chi-square  with  respect  to  W,  evaluated 
at  Wmin((^R),  is  positive  (true  minimum  and  not  a  saddle  or  in¬ 
flection  point). 

The  resulting  (Wmin((pR),  <^r)  ambiguities  are  then  ranked 
by  chi-square,  with  the  “first  rank”  solution  having  the  lowest 
chi-square.  It  should  be  noted  that  the  U%  through  U4  regressions 
have  significantly  larger  retrieval  errors  than  the  first-stage  W 
regression,  and  therefore  the  Wm\u  ’ s  are  inferior  to  the  first- 
stage  W  regression  result.  We  therefore  use  the  first-stage  W 
retrieval  as  the  final  retrieval  result. 

IV.  Median  Filter 

Due  to  regression  errors,  the  first  rank  (most  likely)  solution 
is  not  always  the  closest  ambiguity  to  the  true  wind  direction. 
Improvement  in  ambiguity  selection  beyond  the  first  rank  result 
is  made  in  retrieval  postprocessing,  where  we  apply  a  vector  me¬ 
dian  filter  (MF)  [10]  to  correct  isolated  ambiguity  selection  er¬ 
rors  while  retaining  the  large  scale  discontinuities  (for  example, 
fronts)  found  in  natural  wind  fields. 

The  MF  is  an  iterative  spatial  filter  that  selects  from  the  set 
of  ambiguities  in  the  cell  being  filtered  the  one  most  consistent 
with  the  previously  selected  ambiguities  in  a  window  of  sur¬ 
rounding  cells.  It  accomplishes  this  by  minimizing  the  vector 
median  filter  function  [10],  a  weighted  sum  of  norms  of  vector 
differences  between  an  ambiguity  and  the  previously  selected 
ambiguities  for  all  cells  in  the  window.  After  the  MF  is  applied 
to  all  ocean  cells  in  an  orbit,  the  newly  selected  wind  field  re¬ 
places  the  previous  wind  field.  The  process  is  iterated  until  the 
wind  field  converges. 

In  our  application,  the  window  is  7  x  7  cells  (in  scan-based 
coordinates)  centered  on  and  including  the  cell  being  filtered. 
Missing  cells  in  the  window  (due  to  scan  edges,  land,  ice,  etc.), 
are  “replaced”  with  the  ambiguities  from  the  center  cell.  Be¬ 
cause  wind  direction  retrieval  performance  improves  with  in¬ 
creasing  W  (Section  VI),  the  term  for  each  cell  in  the  filter 
function  is  weighted  linearly  by  the  retrieved  wind  speed  for 
that  cell;  weight  =  min (0.1  *  W,  1.0)  for  high  confidence  re¬ 
trievals,  and  half  that  for  low-confidence  retrievals  (Section  V) 
and  for  replacements. 

The  MF  is  initialized  with  the  first  rank  wind  vector  from 
each  ocean  retrieval  cell.  Optionally,  it  can  be  initialized  with 
nudged  wind  vector  fields  [1 1];  the  nudged  field  is  obtained  by 
selecting  the  first  or  second  rank  ambiguity  for  each  cell  that  is 
closest  to  the  wind  vector  obtained  from  the  NCEP  Global  Data 
Assimilation  System  (GDAS)  1°  x  1°  analysis  closest  in  time, 
spatially  interpolated  to  the  cell  locations. 

V.  Algorithm  Training  and  Testing 

Six  months  of  WindSat  version  1.6.1  Sensor  Data  Records 
(SDRs),  spanning  a  period  when  thermal  gradients  in  the  warm 
load  [12]  are  minimal  (September  2003  through  February  2004), 
have  been  collocated  with  various  sources  of  NWP  and  satellite 
retrievals  for  use  in  regression  training  and  retrieval  testing  (al¬ 
ternating  two  days  training,  one  day  testing  throughout). 

The  training  dataset  consists  of  collocations  between  Windsat 
SDRs,  the  QuikSCAT  science  product  [1 1]  for  wind  vectors,  and 


SSM/I.  GDAS  1°  x  1°  analyses  spatially  interpolated  to  the  lo¬ 
cation  of  Windsat  measurements  are  used  for  SST,  and  for  wind 
vectors  when  a  QuikSCAT  observation  is  not  available.  The  col¬ 
location  with  SSM/I  is  included  to  facilitate  development  of  V 
and  L  regressions.  The  wind  vector  testing  dataset  is  a  two-way 
collocation  between  WindSat  and  QuikSCAT.  For  both  training 
and  testing,  the  spatial  collocation  window  is  25  km  for  all  data 
sources,  and  the  temporal  collocation  window  is  60  min  for 
GDAS  and  QuikSCAT  and  35  min  for  SSM/I.  Before  colloca¬ 
tion,  the  eight  retrieval  cells  along  both  edges  of  the  QuikSCAT 
swath  were  removed,  because  they  contain  less  than  the  optimal 
four  beam  combinations,  and  have  degraded  wind  vectors  [11]. 

The  regression  training  dataset  is  quality  controlled  through 
exclusion  of  every  condition  listed  in  the  next  paragraph  and 
SDRs  that  are  likely  influenced  by  thermal  gradients  in  the  warm 
load.  No  such  exclusions  are  made  in  the  testing  dataset;  instead, 
quality  control  flags  are  set  during  retrieval. 

Retrievals  are  performed  for  all  SDRs  except  those  in  the  aft 
scan  [1],  with  surface  type  other  than  ocean  (no  coast  or  near 
coast),  with  one  or  more  Tg  outside  the  physical  bounds  for 
no  or  light  rain,  or  with  any  EIA  outside  the  expected  range 
(i.e.,  more  than  0.5°  from  nominal).  The  following  are  flagged 
as  “low  confidence  retrievals”:  1)  in  lakes  or  inland  seas;  2)  re¬ 
trieval  cell  may  contain  rain;  3)  in  areas  found  to  periodically 
suffer  from  radio-frequency  interference  (RFI)  or  ice  contam¬ 
ination;  and  4)  are  likely  to  suffer  land  contamination  (within 
105  km  or  1.75  retrieval  cell  HPBW  from  nearest  land).  Condi¬ 
tion  (2)  is  determined  using  an  S  SMI-type  rain  flagging  adapted 
to  WindSat  frequencies  [13].  Static,  geographical  masks  define 
the  other  conditions. 

VI.  Retrieval  Performance 

The  retrieval  performance  in  this  section  excludes  low  confi¬ 
dence  retrievals,  which  are  20%  of  the  total.1  The  Faraday  ro¬ 
tation  correction  necessary  to  optimize  ip  retrieval  performance 
is  under  validation,  and  has  not  been  applied  to  the  SDRs  in  the 
regression  testing  and  training  datasets. 

Wind  vector  retrieval  errors  are  computed  relative  to  collo¬ 
cated  QuikSCAT  data  (Section  V).  The  QuikSCAT  testing  data 
has  errors  relative  to  in  situ  ground  truth,  as  does  the  QuikSCAT 
and  GDAS  data  used  for  wind  vector  regression  training.  How¬ 
ever,  the  errors  computed  here  should  be  good  approximations 
(and  likely  upper  bounds)  to  those  that  would  be  obtained  by 
training  and  testing  with  actual  ground  truth  data  (or  by  training 
with  Tg  simulated  using  a  physical  model  calibrated  to  ground 
truth,  and  testing  with  ground  truth). 

Table  II  shows  the  overall  regression  retrieval  errors.  Given 
the  considerations  above,  the  overall  W  rms  error  is  well  within 
1.0  m/s.  Training  the  W  regression  on  a  natural,  nonuniform 
distribution  of  wind  speeds  (as  contained  in  the  collocated  data) 
initially  resulted  in  a  large  W  bias  error  at  high  wind  speeds 
(—4  m/s  at  23  m/s).  We  greatly  reduced  this  bias  (at  the  expense 

x89%  of  low  confidence  retrievals  were  flagged  for  possible  rain,  3%  for 
RFI/ice,  10%  for  land  contamination,  and  1%  for  inland  lakes  and  seas;  mul¬ 
tiple  flags  can  be  set  for  each  retrieval.  The  rain  flag  was  developed  for  filtering 
all  rain  cases  from  the  training  data,  and  therefore  errs  on  the  cautious  side  and 
overestimates  the  occurrence  of  rain. 
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TABLE  II 

Overall  Retrieval  Errors  for  Regressions 


QuikSCAT  Wind  Speed  (m/s) 


Fig.  1.  WindSat  wind  speed  errors  stratified  by  QuikSCAT  wind  speed. 


of  introducing  a  small  overall  bias)  by  weighting  each  of  the  col¬ 
located  observations  in  the  regression  training  by  the  reciprocal 
of  the  square  root  of  the  wind  speed  distribution  [14].  The  strat¬ 
ification  of  the  W  retrieval  errors  versus  the  QuikSCAT  wind 
speed  is  shown  in  Fig.  1. 

The  overall  regression  errors  for  U\  through  U4  in  Table  II 
are  much  larger  than  those  for  W.  This  occurs  because  the  wind 
speed  signal  in  Th  is  large  [5],  while  the  directional  signals 
are  relatively  small  [6],  [15].  The  U\  through  U4  regressions 
must  extract  parts  of  the  small  directional  signals,  while  com¬ 
pensating  for  atmospheric  and  SST  variations.  Because  of  the 
sine/cosine  dependence  of  the  wind  direction  signals,  the  U\ 
and  C/3  regressions  use  mainly  Tv  and  Th,  while  the  U2  and  C/4 
regressions  use  mainly  Tjj  and  T4.  Because  the  wind  direction 
signals  for  Tv  and  Th  are  superposed  on  a  large  and  highly  vari¬ 
able  nondirectional  component,  the  regression  errors  for  U\  and 
C/3  are  larger  than  those  for  C/2  and  C/4. 

Fig.  2  shows  the  wind  direction  rms  errors  for  the  first  ranked 
(FR),  median  filtered  (MF),  median  filtered  with  nudging 
(MFNG)  ambiguities,  and  the  closest  ambiguity  to  the  collo¬ 
cated  wind  direction  (CL).  The  wind  direction  bias  error  at  all 
W  is  small  (less  than  3°)  and  is  not  shown. 

Above  12  m/s,  the  agreement  between  the  wind  direction  re¬ 
trievals  and  QuikSCAT  is  excellent  for  all  of  the  above  ambi¬ 
guities.  Our  initial  wind  direction  performance  is  within  20°  of 
QuikSCAT  for  wind  speeds  above  9.9,  8.8,  7.6,  and  7.2  m/s  for 
the  FR,  MF,  MFNG,  and  CL  ambiguities,  respectively.  Below 
5  m/s,  the  wind  direction  signals  are  very  small  [6],  [15],  and 
the  wind  direction  retrieval  performance  degrades  accordingly. 

Fig.  3  shows  the  skill  (percentage  of  retrievals  where  the  se¬ 
lected  ambiguity  is  the  CL  ambiguity)  for  the  FR,  MF,  and 
MFNG  ambiguities,  and  the  total  number  of  ambiguities.  De¬ 
spite  high  FR  skill,  the  median  filter  by  itself  makes  only  a 
small  improvement  in  the  wind  direction  error  (6°  to  8°  below 
8  m/s,  less  than  3°  above).  One  reason  is  that  the  number  of 


Fig.  2.  WindSat  wind  direction  rms  errors  stratified  by  QuikSCAT  wind  speed. 
Errors  are  determined  relative  to  the  collocated  QuikSCAT  science  product. 


Fig.  3.  Wind  direction  skill  and  mean  number  of  ambiguities  stratified  by 
QuikSCAT  wind  speed. 


ambiguities  is  low.  Due  to  the  requirement  that  the  MLE-re- 
trieved  wind  speed  be  nonnegative  (Section  III),  there  are  no 
more  than  two  ambiguities  for  any  retrieval,2  with  an  average  of 
1.56  over  all  retrievals.  Thus,  the  median  filter  cannot  improve 
44%of  the  retrievals.  Another  reason  is  that,  when  one  considers 
only  retrievals  with  two  ambiguities,  the  skill  is  reduced  from 
that  shown  in  the  figure;  for  example  the  FR  skill  drops  from 
84%  to  76%  at  5  m/s. 

Additionally,  MF  wind  fields  show  good  qualitative  agree¬ 
ment  with  QuikSCAT  wind  fields  except  inside  some  regions 
where  W  <  6  m/s  and  wind  direction  is  changing  rapidly;  there, 
low-skill  degraded  retrievals  often  occur  in  contiguous  patches. 
By  itself,  the  median  filter  rarely  improves  upon  the  FR  am¬ 
biguity  in  such  patches.  Analogous  to  scatterometer  retrievals, 
we  have  found  that  the  separation  between  the  first  and  second 
rank  ambiguities  is  about  180°,  which  may  further  inhibit  me¬ 
dian  filter  performance  in  and  near  the  patches.  Nudging,  which 
can  select  the  second  rank  ambiguity  to  initialize  the  median 

2Without  this  requirement,  the  maximum  number  of  ambiguities  is  3.  How¬ 
ever,  ambiguities  with  negative  IUmin  always  have  the  highest  chi-square  of 
any  ambiguity  (often  by  an  order  of  magnitude),  have  a  value  near  —1  m/s, 
regardless  of  the  magnitude  of  the  wind  speeds  computed  from  (Uly  U2)  and 
(U3,  U4)  vectors,  and  appear  when  these  two  wind  speeds  are  grossly  incon¬ 
sistent.  Thus,  we  consider  ambiguities  with  negative  Wm in  to  be  an  unphysical 
result  obtained  when  the  MLE  is  given  inconsistent  values  of  U±  through  U 4  . 
As  such,  they  should  be  discarded. 
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filter,  allows  significant  improvement  in  the  agreement  between 
MFNG  and  QuikSCAT  wind  fields  over  the  patches,  and  thus 
in  wind  direction  performance  at  low  and  intermediate  wind 
speeds. 

Because  nudging  chooses  between  two  ambiguities  separated 
by  approximately  180°,  small  to  moderate  perturbations  to  the 
wind  field  used  for  nudging  will  have  little  effect  on  the  nudged 
wind  field  used  to  initialize  the  median  filter  [11].  Therefore,  the 
incorporation  of  QuikSCAT  data  into  GDAS  wind  fields  should 
not  significantly  affect  the  nudged  wind  field,  or  the  comparison 
of  the  MFNG  ambiguities  with  QuikSCAT. 

VII.  Conclusion 

We  have  developed  a  regression  retrieval  technique  for 
WindS  at  retrieval  of  ocean  surface  wind  speed  and  a  com¬ 
bination  regression/MLE  technique  for  ocean  surface  wind 
direction  retrieval.  The  regressions  were  derived  empirically, 
using  collocated  NWP  data  and  satellite-based  retrievals.  To¬ 
gether,  they  provide  an  initial  assessment  of  WindSat  retrieval 
capabilities.  Comparisons  of  the  initial  WindSat  retrievals  with 
QuikSCAT  wind  vectors  have  demonstrated  the  capability  for 
passive  microwave  wind  direction  retrieval,  and  show  excellent 
agreement  with  QuikSCAT  wind  speeds. 

We  expect  the  retrievals  to  improve  due  to  ongoing  work  on 
the  SDR  algorithm,  quality  control  flagging,  and  the  retrieval 
algorithms.  Ongoing  work  on  the  SDR  algorithm  includes  ad¬ 
dition  of  corrections  for  Faraday  rotation  and  thermal  gradi¬ 
ents  in  the  warm  calibration  load.  We  are  also  working  on  im¬ 
proved  rain  flagging,  and  replacement  of  the  sea  ice  and  RFI 
geographical  masks  with  dynamic  ice  and  RFI  detection.  Re¬ 
trieval  algorithm  work  includes  using  nonuniform  bin  widths  in 
the  second-stage  regressions,  tuning  the  median  filter,  and  in¬ 
corporating  buoy  observations  into  the  regression  training  and 
testing. 
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